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ABSTRACT 

We report the discovery of a planet transiting a moderately faint (V = 12.3 mag) late F star, 
with an orbital period of 3.92289 ± 0.00004 days. From the transit light curve and radial velocity 
measurements we determine that the radius of the planet is R p = 1.40 ± 0.06i?j up and that the 
mass is M p = 0.78 ± 0.09 Mj up . The density of the new planet, p p — 0.35 ± 0.06 gem -3 , fits to the 
low-density tail of the currently known transiting planets. We find that the center of transit is at 
T c = 2,454,417.9077 ± 0.0003 (HJD), and the total transit duration is 0.143 ± 0.004days. The host 
star has M* = 1.28 ± 0.13 Mq and i?* = 1.32 ± 0.07 R Q . 

Subject headings: stars: individual: GSC 02463-00281 - planetary systems: individual: HAT-P-9b 



1. INTRODUCTION 

Transiting extra-solar planets are important astrophys- 
ical objects as they allo w to test planetary structure 
and evolution theory (e.g . jFortn cv 2008a; Bu rrows et al.l 
120081: iBaraffe et alfl2008l) especially because they yield 
a measurement of the planetary mass and radius. 



Over the last few years, the sample of known tran- 
siting planets has grown substantially, leading to an 
improved th eoretical understanding of their physical 
nature (e.g. 



iFortnev et al.1 



Guillot et all 120061: iBurrows et al. 1 120071 : 



i hoOSbHChabrier fc Baraffell2007l L In addi- 
tion, the increasing sample of transiting planets enabled 
the discovery of some int eresting correlations, such as 
the mass-period r elation (|Mazeh. Zucker. fc Pontl [20051 : 
iGaudi et al.ll2005f l. The astrophysics behind these corre- 
lations is not fully understood, and a clear way towards 
progress is the discovery of many more transiting planets. 

We report here the discovery of anoth er transiting 
planet dete cted by the HATNet project 11 (iBakos et al.l 
120021 12004D . labeled HAT-P-9b, and our determination 
of its parameters, including mass, radius, density and 
surface gravity. In § [2] we describe our photometric and 
spectroscopic observations and in § [3] we derive the stel- 
lar parameters. The orbital solution is performed in § [4] 

1 Based in part on radial velocities obtained with the SOPHIE 
spectrograph mounted on the 1.93m telescope at the Observatory 
of Haute Provance (runs 07A.PNP.MAZE, 07B.PNP.MAZE, 
08A.PNP.MAZE). 
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11 http://www.hatnet.hu 




••• • •• J> • - v , 

■s 099 ■ t *. . ' 



98 

:; ; ;? 



Fig. 1. — The top panel shows the unbinned HATNet light 
curve with 6884 data points, phased with the orbital period of 
P = 3.92289 d. The binned light curve is over-plotted. The center 
of the transit is at phase zero. The bottom panel shows a zoom- 
in around the transit. The HATNet light curve is folded with 
the e phem eris obtained here, based on the follow-up light curves 
(see § 12. 4jl . The small offset between phase zero and transit center 
suggests there is a slight difference in the ephemeris for the HATNet 
data and the photometric follow-up data, obtained ~ 1000 days 
apart. See §[7]for further discussion. 

and a discussion of possible blend scenarios is brought 
in § O The determination of the light curve parameters 
and the planet's physical parameters is described in §[5] 
and we bring a discussion in § [7j 

2. OBSERVATIONS AND ANALYSIS 

2.1. Detection of the transit in the HATNet data 

HAT-P-9 is positioned in HATNet 's internally labeled 
field G176, centered at a = 07 h 28 h , 5 = 37°30'. This 
field was observed in network mode by the HAT-6 tele- 
scope, located at the Fred Lawrence Whipple Observa- 
tory (FLWO) of the Smithsonian Astrophysical Observa- 
tory (SAO), and the HAT-9 telescope at the Submillime- 
ter Array (SMA) site atop Mauna Kea, Hawaii. In total, 
6884 exposures were obtained at a 5.5 minute cadence 
between 2004 November 26 and 2005 October 21 (UT). 

Preliminary reduction included standard bias, dark 
and flat-field corrections, followed by an astrometric solu- 

Photometry 



tion using the code of iPal fc Bako; 



>v an asti 
i (120061) . 
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was applied using fine-tuned aperture photometry while 
the raw light curves were processed by our external pa- 
rameter decorrelatio n technique. Next, the Trend F ilter- 
ing Algorithm (TFA; IKovacs. Bakos. fc Novesll2005h was 
applied to get the final light curves. To search for the sig- 
nature of a transiting planet i n the light curves we used 
the B ox Least Squares (BLS; IKovacs. Zucker. fc Mazehl 
I2002f) algorithm. 

A transit-like signal was identified in the light 
curve of HAT-P-9 (GSC 02463-00281, 2MASS 
J07204044+3708263), which is a U=12.3mag star, 
fainter than most transiting planet host stars detected 
with small-aperture, wide-field, ground-based cam- 
paigns. The HATNet light curve is presented at the top 
panel of Fig. [1] folded on the ephemeris obtained from 
analyzing the follow-up light curves (P — 3.92289 days 
and T c = 2454417.9077, see 6ETlj). showing a flux 
decrement of ~ 1 % at phase zero. The figure shows 
a small shift between phase zero and transit center, 
suggesting a possible slight difference in the ephemeris 
for the HATNet data and the photometric follow-up 
data, obtained ~ 1000 days apart (see § [7] for further 
discussion). This transiting planet candidate, along with 
others from the same field, was selected for follow-up 
observations, to investigate its nature. 

2.2. Early spectroscopy follow-up 

Initial follow-up observ ations were m ade with the CfA 
Digital Speedometer (DS; [Latham 1993) in order to char- 
acterize the host star and to reject obvious astrophysi- 
cal false-positive scenarios that mimic planetary tran- 
sits. The five radial velocity (RV) measurements ob- 
tained over an interval of 62 days showed an rms resid- 
ual of 1.1 km s -1 , consistent with no detectable RV vari- 
ation. Atmospheric parameters for the star (effective 
temperature T c g, surface gravity logg, and projected 
rotational velocity fsini) were derive d as described by 
iTorres. Neuhauser fc Guentherl (|2002f ). initially assum- 
ing a fixed metallicity of [Fe/H] = 0.0 ± 0.2. We ob- 
tained log g = 4.05 ± 0.30, T cff = 6130 ± 150 K and 
vs'mi = 12.2±1.0kms _1 . As the DS results were consis- 
tent with a planet orbiting a moderately-rotating main 
sequence star, this target was selected for high-precision 
spectroscopy follow-up. 

2.3. High-precision spectroscopy follow-up 

Observations were carried out at the Haute Provence 
Observatory ( OHP) 1.93-m telescope, with the S OPHIE 
spectrograph (jBouchv fc the Sophie Team! 120061) . SO- 
PHIE is a multi-order echelle spectrograph fed through 
two fibers, one of which is used for starlight and the other 
for sky background or a wavelength calibration lamp. 
The instrument is entirely computer-controlled and a 
standard data reduction pipeline automatically processes 
the data upon CCD readout. RVs are calculated by nu- 
merical cross-correlation with a high resolution observed 
spectral template of a G2 star. 

HAT-P-9 was observed with SOPHIE in the high- 
efficiency mode (R ~ 39000) during three observing 
runs, from 2007 May until 2008 May. Due to the rel- 
ative faintness of this star, exposure times were in the 
range of 25 to 75 minutes, depending on observing con- 
ditions. The resulting signal to noise ratios were 32-46 
per pixel at A = 5500 A. Using the empirical relation 
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Fig. 2. — The top panel shows the RV measurements phased 
with the period of P = 3.92289 days and mid-transit time of 
T c = 2454417.9077 (HJD). The zero-point in phase corresponds 
to the epoch of mid-transit. Overlaid is the best sinusoidal fit. 
The bottom panel shows the residuals from the fit. 



TABLE 1 
Radial Velocities for HAT-P-9. 



BJD -2,400,000 


RV a 


"FtV 


BS b 


S/N c 


(days) 


(ms _1 ) 


(ms- 1 ) 


(ms- 1 ) 




54233.3347 


22667 


24 


-62 


43 


54379.6357 


22586 


20 


5 


14 


54380.6213 


22672 


24 


-93 


38 


54381.5815 


22749 


23 


-148 


40 


54413.6292 


22727 


21 


3 


43 


54414.6084 


22610 


24 


-31 


38 


54415.6457 


22675 


27 


-45 


34 


54416.6392 


22740 


21 


-13 


42 


54421.6701 


22659 


31 


-11 


32 


54422.6651 


22595 


20 


-59 


46 


54587.3543 


22536 


21 


10 


14 


54588.3571 


22650 


23 


30 


41 


54589.3525 


22730 


23 


-2 


40 


54590.3425 


22699 


25 


-45 


38 


54591.3647 


22564 


30 


10 


32 



a The RVs include the barycentric correction. 

Bisector span. 
c Signal to noise ratio per pixel at A — 5500 A. 



of iCameron et all (|2007f ) we estimated the RV photon- 
noise uncertainties to be 20-31 ms" 1 . We made 17 RV 
measurements in total. Two of those were highly con- 
taminated by the Moon and were ignored, leaving 15 RV 
measurements, listed in Table [T] 

2.4. Photometry follow-up 

In order to better characterize the transit parameters 
and to derive a better ephemeris, we performed photo- 
metric follow-up observations with 1-m class telescopes. 
We obtained a total of four transit light curves of HAT- 
P-9b, shown in Fig. [3] Two events were observed by the 
KeplerCam detecto r on the FLWO 1.2 m telescope (see 
iHolman et a!1l2007t ) on UT 2007 Nov 13 and UT 2008 
Jan 11, in the Sloan z-band. We refer to the 2007 Nov 
13 event as having a transit number N tI — 0, so the 2008 
Jan 11 transit number is N tr — 15. In addition, two 
light curves were obtained at the Wise Observatory. On 
UT 2007 Dec 6 the AT tr = 6 t ransit event was ob served 
by the Wise 0.46 m telescope (jBrosch et al J 12008). with 
no filter. The following event, with N tI — 7, was ob- 
served on UT 2007 Dec 10 by the Wise 1 m telescope in 
the Cousins i?-band. An additional light curve, obtained 
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TABLE 2 

List of follow-up light curves of HAT-P-9. 





Start Date 


Observatory+ 


Filter 




U2 


Cadence 


/3 a 


RMS 




UT 


Telescope 








min - 1 




% 





2007 Nov 13 


FLWO 1.2 m 


z 


0.1313 


0.3664 


1.0 


1.1 


0.17 


6 


2007 Dec 6 


Wise 0.46 m 


clear 


0.2403 


0.3816 


1.1 


1.0 


0.21 


7 


2007 Dec 10 


Wise 1.0 m 


R 


0.2403 


0.3816 


0.1 


1.2 


0.22 


15 


2008 Jan 11 


FLWO 1.2 m 


z 


0.1313 


0.3664 


1.0 


1.6 


0.18 



The correlated noise factor, by which the errors of each light curve are multiplied (see § 12. 4t . 




N =7 Wise 1m ft 



N,=15 FLWO 1.2m z 



TABLE 3 

Summary of stellar parameters for HAT-P-9. 



-0.1 -0.05 0.05 0.1 0.15 

HJD-T [days] 



Fig. 3.— Follow-up light curves of HAT-P-9. The four light 
curves are shown at the top of the figure, with the fitted model 
over-plotted. The label next to each light curve gives the transit 
event number, observatory, telescope and filter used. Residuals are 
presented at the bottom of the figure, with the same top to bottom 
order as the actual light curves. 



with the FLWO 1.2 m telescope, was of poor quality and 
is not included here. Table [2] lists for each light curve its 
Nn number, UT date, observatory, telescope, filter used, 
limb-darkening coefficients used in its analysis, mean ca- 
dence, the correlated noise (3 factor (see below) and the 
RMS residuals from the fitted light curve model. 

Data were reduced in a similar manner to the HAT- 
Net data, using aperture photometry and an ensem- 
ble of ~100 comparison stars in the field. An analytic 
model was fitted to these data, as described below in § 
and yielded a period of 3.92289 ± 0.00004 d and a refer- 
ence epoch of mid-transit T c = 2,454,417.9077 ± 0.0003d 
(HJD). The length of the transit as determined from this 
joint fit is 0.143±0.004d (3 hours, 26 minutes), the length 
of ingress is 0.019 ± 0.003d (27 minutes), and the central 
transit depth is 1.17 ± 0.01 %. The latter value is simply 
the square of the radius ratio (see Table [4|) if one ignores 
the limb-darkening effect. This effect increases the depth 
by about 0.1 %, in the bands used here. 



Parameter 


Value 




Source 




R.A. 


07 h 20 m 40? 44 


2MASS 








Dec. 


+37° 08' 26/3 


2MASS 








m v (mag) 


12.297 ±0.063 


TASS 








T cff (K) 


6350 ± 150 


DS ± Yonsei-Yale 






v s'mi (km s ) 


11.9 ± 1.0 


DS + Yonsci-Yale 






logs 


4 OQ+ 03 


DS + Yonsei-Yale + 


light 


curve shape 


[Fc/H] (dex) 


0.12 ±0.20 


DS + Yonsei-Yale 






Mass (M ) 


1.28 ±0.13 


Yonsei-Yale - 


h light 


curve 


shape 


Radius (Rq) 


1.32 ±0.07 


Yonsei-Yale - 


h light 


curve 


shape 


log(WL Q ) 


41 + ' 08 


Yonsei-Yale 








My (mag) 


, 7 +0.3 

— 0.2 


Yonsei-Yale 








Age (Gyr) 


1 6+ 18 

t.D_l.4 


Yonsei-Yale - 


h light 


curve 


shape 


B-V (mag) 


<50+°- 0a 


Yonsei-Yale 








Distance (pc) a 


480 ± 60 


Yonsei-Yale 









a Assuming extinction of A{V) — 0.15 mag, see text at § [3l 

3. STELLAR PARAMETERS 

The mass (M p ) and radius (Rp) of a transiting planet, 
determined from transit photometry and RV data, is 
dependent also on those of the parent star. In or- 
der to determine the stellar properties needed to place 
M p and R p on an absolute scale, we made use of stel- 
lar evolution models along with the observational con- 
st raints from spectros copy and photometry, as described 
in lTorres et al.l (|2008f) . Because of its relative faintness, 
the host star does not have a parallax measurement from 
Hipparcos, and thus a direct estimate of the absolute 
magnitude is not available for use as a constraint. An 
alternative approach is to use the surface gravity of the 
star, which is a measure of the evolutionary state of the 
star and therefore has a very strong influence on the ra- 
dius. However, log g is a difficult quantity to measure 
spectroscopically and is often strongly correlated with 
other spectros copic parameters (see § 12. 2p . It has been 
pointed out bv lSozzetti et al.1 (|2007l ) that the normalized 
separation of the planet, a/i?*, can provide a much bet- 
ter constraint for stellar parameter determination than 
the spectroscopic log g. The a/R* quantity can be deter- 
mined directly from the photometric observations with 
no additional assumptions (other than the values of the 
parameters describing limb-darkening, which is a second- 
order effect), and it is related to the mean density of the 
host star. As discussed later in § [6l an analytic fit to 
the light curve yields a/i?+ = 8.6 ±0.2. This value, along 
with T c g from § 12. 2\ and assuming [Fe/H] = 0.0±0.2, was 
comp ared with the Y onsei-Yale stellar evolution mod- 
els of lYi et al.l (|2001f ) . This resulted in values for the 



stellar mass and radius of M*= 1.16 
1.28 ±0.08 i?©, and an estimated age of 3.6 



+0.12 
-0.15 



M Q and R*= 
+^Gyr. The 
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result for logo was 4.29^q'q^, consistent with the value 
derived from the DS spectra. 

Nevertheless, in order to verify the overall consistency 
of our results and refine the stellar parameters, we carried 
out a new iteration. We imposed this latter value of log g 
(coming from stellar evolution modeling), and analyzed 
the DS spectra by allowing the metallicity, vs'mi and 
T off to vary. The new results were [Fe/H] = 0.12 ± 0.20, 
usini = 11.9 ± 1.0 km s -1 and T cS = 6350 ± 150 K. Re- 
peating the stellar evolution modeling resulted in M*= 
1.28 ± 0.13 M©, R k = 1.32 ± 0.07 R e and an age of 
1.6l^?Gyr. The stellar properties are summarized in 
Table [3] and they correspond to a late F star. 

We used also our SOPHIE high-resolution spectra 
to estimate usini and [Fe/H]. Based on our result 
for B — V from the stellar evolution models, of B — 
V = 0.50tom mag, we got vsini = 10.1 ± l.Okms" 1 
and [Fe/H]= +0.18 ± O.lOdex. Those values are close 
to the results from the DS spectra and the stellar evolu- 
tion model. However, the SOPHIE spectra were taken in 
High-Efficiency mode, where there is a known problem 
with removing the echelle blaze function. Hence, these 
estimates are used only for comparison, and are not in- 
cluded in our final result. 

To check our spectroscopically determined T c g we 
used several publicly avai lable color indices for this 
star, and the calibration s of iRarmrez fe Melendezl ([2005h 
and ICasagrande et al.l (|2006h to derive independent 
temperature estimates, using the [Fe/H] value found 
above in these calibrations. We adjusted the redden- 
ing till the photometric temperature matched the spec- 
troscopic temperature, yie lding E(B — V) = 0.053 
([Ramirez fc Melendezl [20 05) and E(B - V) = 0.043 
(|Casagrande et aLll2006f) . We adopted the mean value 
of 0.048 as reddening, im plying an extinction of A( V) = 
0.15 mag. Note that the iBurstein fc Heilesl (| 19821 ) red- 
dening maps for this celestial posit ion (/, b = 181. 1 , 21.4) 
give E(B - V) = 0.072 and the ISchleeel et all (fl998h 
maps yield E(B — V) = 0.065, broadly consistent with 
our findings, especially if we take into account that these 
latter methods measure the total reddening along the line 
of sight. 

Using the V m agnitude from the TASS survey 
(jDroege et al.ll2006D . V = 12.297 ± 0.063 mag, the abso- 
lute V magnitude from the stellar evolution model, and 
the A(V) — 0.15 mag value determined above, the dis- 
tance to HAT-P-9 is 480 ± 60 pc. 

4. SPECTROSCOPIC ORBITAL SOLUTION 

Our 15 RV measurements from SOPHIE were fitted 
with a Keplerian orbit model solving for the velocity 
semi-amplitude K and the center-of-mass velocity 7, 
holding the period and transit epoch fixed at the well- 
determined values from photometry (see Table [4]). The 
eccentricity was initially set to zero. The fit yields 
K = 84.7 ± 7.9 ms" 1 and 7 = 22.665 ± 0.006 km s" 1 . 
The observations and fitted RV curve are displayed in 
the top panel of Fig. O The residuals are presented in 
the bottom panel of the same figure. RMS residuals is 
22.1ms _1 , consistent with the RV uncertainties. The 
value of x 2 is 12.4 for 13 degrees of freedom. 

5. EXCLUDING BLEND SCENARIOS 
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Fig. 4. — Line bisector spans folded on the orbital phase. Albeit 
one or two outliers, the bisector spans do not show a variation with 
the orbital phase. 

We tested the reality of the velocity variations by ex- 
amining the spectral line bisector spans (BSs) of the 
star using our SOPHIE data. If the measured velocity 
changes are due only to distortions in the line profiles 
arising from contamination of the spectrum by the pres- 
ence of a binary with a period of 3.92 days, we would 
expect the BSs (which measure line asymmetry) to vary 
with this period, resulting in a correlation between B S 
and RV (see, e.g.. lQueloz etHll2"o5H : iTorres et al.ll2005l) . 
As shown in Fig. |U the BSs show no significant vari- 
ations, except one or two outliers. The correlation co- 
efficient between BS and RV for all 15 measurements is 
-0.38. Ignoring the extreme point reduces the correlation 
to -0.18. 

In order to estimate the statistical significance of a 
correlation between BS and RV we define the following 
statistics: 

r a = (1) 

CTBS 

where ass is the standard deviation of the BS values 
and gbs, fit the standard deviation of the residuals of 
a BS fit to the RVs. Both standard deviations are the 
unbiased estimators. For pure noise r a will equal 1.0. 
When fitting the BSs to the RVs of all 15 measurements, 
using polynomials of degrees 1-3 we got r a in the range 
0.984-1.043, consistent with no significant correlation. 

Another sign of a binary would be a dependence be- 
tween the RV amplitude and the template used. This 
may happen when the components of the blended binary 
are of a different spectral type than the primary target. 
We re-calculated the RVs using F0 and K5 templates 
and got an amplitude consistent with our original value, 
within lcr (not shown). 

These analyses indicate that the orbiting body is a 
planet with high significance. 

6. PLANETARY PARAMETERS 

To determine the light curve parameters of HAT-P-9b 
we fitted the four availa ble transit ligh t curv es simulta- 
neously, as described in iShporer et a l. (2008). A circu- 
lar orbit was assumed. We adopted a quadratic limb- 
darkening law for the star, and determined the ap propri- 
ate coefficien ts, u\ and U2, by interpolating on the lClaretl 
(|2000L l2004h grids for the atmospheric model described 
above. For the Wise 0.46 m light curve, observed with no 
filter, we adopted the Cousins i?-band limb-darkening co- 
efficients as the CCD response resembles a "wide-i?" filter 
(jBrosch et al.l [2008^1 . The values of the limb-darkening 
coefficients used for each light curve are listed on Table [H 

The drop in flu x in the light curves wa s modeled with 
the formalism of iMandel fc Agol (|2002f ). The five ad- 
justed parameters in the fit were i) the period P\ ii) 
the mid-transit time of the N tr = transit event, T cU , 
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TABLE 4 

Orbital fit and planetary parameters for 
the hat-p-9b system. 



Parameter 


Value 


Ephcmeris: 




Period (day) a 


3.92289 ± 0.00004 


T c (HJD) a 


2,454,417.9077 ± 0.0003 


Transit duration (day) 


0.143 ±0.004 


Ingress duration (day) 


0.019 ±0.003 


Orbital parameters: 




7 (kms -1 ) 


22.665 ± 0.006 


K (ms" 1 ) 
e a 


84.7 ± 7.9 



Light curve parameters: 




a/R* 


8.6 ±0.2 


Rp/ R+ 
b 


0.1083 ±0.0005 


0.52 ±0.03 



i (deg) 
a (AU) 
M p (Mj) 
Rp{Rj)_ 
' ; cm 3 ) 

-2\b 



Pp 

9 P (ms" 



2 )' 



d 

f p (10 9 erg cm" 



86.5 ±0.2 
0.053 ± 0.002 
0.78 ±0.09 
1.40 ±0.06 
0.35 ± 0.06 
9.8 ± 1.0 
1530 ± 40 
0.046 ± 0.007 

1 Q + 0.3 



Fixed in the orbital fit. 



quantities, see 



Based on only directly observable 
ISouthworth et al.l J2007I , Eq. 4). 
c Planetary thermal-equilibrium surface temperature. 
^ Safronov number, sec Hansen &: Barman ( 2007. Eq. 
2) 

G Stellar flux at the planet. 

denoted simply as T c ; iii) the relative planetary ra- 
dius, r = Rp/Ri,; iv) the orbital semi-major axis scaled 
by the stellar radius, a/R+; v) the impact parameter, 
b = acos(i)/i?,t, where i is the orbital inclination angle. 

Acc ounting for corr elated noise in the photometric 
data (|Pont et al.l [2006) was done si milarly to the "time- 
averaging" method of ?, see also IShporer et all [2008) . 
After a preliminary analysis we binned the residual light 
curves using bin sizes close to the duration of ingress and 
egress, which is 27 minutes. The presence of correlated 
noise in the data was quantified as the ratio between the 
standard deviation of the binned residual light curve and 
the expected standard deviation assuming pure white 
noise. This ratio is calculated separately for each bin 
size and we defined (3 to be the largest ratio among the 
bin sizes we used. For each light curve we multiplied the 
relative flux errors by (3 and repeated the analysis. The 
j3 value of each light curve is listed on Table [5] 

We used the Markov Chain Mont e Carlo algorithm 
(MCMC, see, e.g. lHolman et al.ll2007t ) to derive the best 
fit parameters and their uncertainties. The MCMC al- 
gorithm results in a distribution of each of the fitted pa- 
rameters. We took the distribution median to be the best 
fit value and the values at 84.13 and 15.87 percentiles to 
be the +ler and -la confidence limits, respectively. 

Our four follow-up light curves are presented in Fig. [3] 
where the fitted model is overplotted and the residu- 
als are also shown. The result for the radius ratio is 
Rp/R* = 0.1083±0.0005, and the normalized separation 
is a/i?* = 8.6 ±0.2. From the mass of the star (Table [3|), 
the orbital parameters (§ |4]) and the light curve param- 
eters, the physical planetary parameters (such as mass, 
radius) are calculated in a straightforward way, and are 
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Fig. 5. — Top panel: The radius-mass diagram for the known 
transiting planets, in log-log scale, with HAT-P-9b marked by 
an open circle. Constant density lines are over-plotted. The 
diagram visually shows that the new planet has a relatively low 
density. The Neptune-mass planet GJ 436b, orbiting an M star, is 
positioned beyond the lower-left corner of this diagram. 
Bottom panel: The mass-period diagram for the known transiting 
planets, in linear scale, with HAT-P-9b marked by an open 
circle. Three planets are positioned outside the boundaries of this 
diagram: The long-period planet HD 17156b, and the massive 
planets HAT-P-2b (a.k.a. HD 147506b) and XO-3b. The planet 
with the lowest mass in this diagram is GJ 436b, orbiting an M 
star, and the one with the largest mass is CoRoT-Exo-2b. 
This figure is based on data taken from 
http://www.inscience.ch/transits/ on 2008 June 1st. 

summarized in Tabled] We note that a/i?*, as derived 
from the light curve fit, is an important constraint for the 
stellar parameter determination (§[3]), which in turn de- 
fines the limb-darkening coefficients that are used in the 
light curve analytic fit. Thus, after the initial analytic 
fit to the light curves and the stellar parameter determi- 
nation, we performed another iteration of the light curve 
fit. We found that the change in parameters was imper- 
ceptible. 

7. DISCUSSION 

We present here the discovery of a new transiting extra- 
solar planet, HAT-P-9b, with a mass of 0.78 ± 0.09 Mj, 
radius of 1.40 ± 0.06 Rj and orbital period of 3.92289 ± 
0.00004days. The V = 12.3 mag host star is among the 
faintest planet-host stars identified with small-aperture, 
wide-field, gro und-based campaigns . In fact, together 
with WASP-5 ([Anderson et al.ll2008f ) they are currently 
the faintest planet-host stars discovered with such cam- 
paigns, where the transit depth is below 2%. 

To show the main characteristics of the new planet 
relative to the currently known transiting planets we plot 
in Fig. [5] the position of the latter on the radius-mass 
diagram, on the top panel, and the mass-period diagram, 
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on the bottom panel. In both panels HAT-P-9b position 
is marked by an open circle. 

The radius-mass diagram, presented here in log- 
log scale, visually shows that the new planet has 
a low density, of p p = 0.35 ± 0.06gcm~ 3 , similar 
to that of HP 2 09458b (e.g., iKnutson et al 
Mazeh et all |2000D . WASP-lb (ICameron et al 



2007 



2007 



Shporer et alj|2007t ICharbonneau et all 120071) . HAT-P 
lb dBakos et alJl2007t IWinn et alJl2007ti ) and CoRoT- 
Exo-lb (|Barge et al.ll2008fh 

Comparing the measu red radius to the modeled radius 
of lBurrows et al. I (|2007h . for the same stellar luminosity, 
semi-major axis and (core-less) planet mass as the HAT- 
P-9 system, shows the measured value is larger than the 
modeled one by 1-2 a. In this comparison we assumed 
sol ar opacity for t h e plan etary atmosphere. As shown 
by iBurrows et al. I (|2007t ). the transit radius effect to- 
gether with enhanced planetary atmospheric opacities, 
relative to solar opacity, can be responsible for inflating 
the planetary radius. Since atmospheric opacities are 
related to the host star metallicity, and HAT-P-9 metal- 
licity seems to be super-solar, it is likely that these two 
effects acc ount for the radius dif ference, as is the case for 
WASP-lb (jBurrows et al. 1120071 see their Fig. 7). Other 
mechanisms that may explain the increased radius in- 
volve d^wrrw^xd_Jraris2ort and dissipation of kinetic en- 
ergy dShowman fe Guillotll2002f ) and layered convection 
dChabrier fc Baraffdl2007n . 

According to its T cq and Safronov number, listed on 
Ta ble II HAT P-9b is li kely a Class II planet, as defined 
by lHansen fc Barman] (|2007f ). This classification is in 
accordance with its low density, since members of this 
class usually have smaller masses and larger radii than 
those of Class I. 

The stellar incident flux at the planet is f p = 
1. 31q3 10 9 erg cm~ 2 s -1 , which is at the lower end of 
the f p range for pM planets, according to the pL/pM 
planet classification of iFortnev et al.l (|2008bf ). Since it is 
positioned near the transition region in f p between these 
two classes of planets it may be used to further charac- 
terize them. 

In the mass-period diag ram the new pl a net is 
positioned n ear XO-lb <|McCulloughlt7ai1 [20061: 
iHolman et all l200l ). OGLE-TR-182b (|Pont et all 



20081) OGLE-TR-lllb dPont et all l200l IWinn et all 
2007al) and HAT-P-6b (|Noves et alTl2008f i" The 
position of HAT-P-9b, along with that of neighbour- 
ing planets, suggests that the mass-period re lation 
(jMazeh. Zucker. fc Pontl 120051 : iGaudi et alJ l2005h . i.e., 
the decrease of planetary mass with increasing orbital 
period, levels-off at periods > 3 days. 

The ephemeris obtained here is based on the four 
follow-up light curves. We compared it to the mid-transit 
time of the HATNet light curve, obtained some 1000 days 
earlier, by propagating the error on the period. The re- 
sult shows a difference of about 1.5 a between the mid- 
transit times of the follow-up light curves and that of the 
HATNet light curve. Future follow-up light curves will be 
able to study more thoroughly this possible discrepancy. 

Finally, we note that due to the line-of-sight stellar 
rotation velocity of 11.9 ± 1.0 km s -1 and the non-zero 
impact parameter 0.52 ± 0.03, this planet is a good can- 
didate for t he measurement of the Ross iter-Mclaughlin 
effect ^e.g.. IWinn et al.ll2005l ). Eq. 6 of IGaudi fc WTiml 
( 2007) gives an expected amplitude for the effect of more 
than 100 ms -1 , which is larger than the orbital ampli- 
tude. However, since this is a relatively faint star, with 
V = 12.3 mag, such a measurement will be challenging. 
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